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Abstract  Detection  and  tracking  are  normally  considered  separate  processes. 
First  the  signal  processing  system  associated  with  a  sensor  examines  the  signal 
to  determine  whether  to  call  detection.  Once  detection  is  called,  it  is  converted 
to  an  estimate  of  one  or  more  of  the  components  of  the  target's  kinematic  state, 
e.g.,  bearing,  position,  or  velocity.  This  estimate  (contact)  is  sent  to  a  tracking 
system  that  determines  whether  the  contact  should  be  associated  with  an 
existing  track  or  used  to  generate  a  new  one.  This  process  works  well  in  high 
signal-to-noise  ratio  (SNR)  situations  but  sacrifices  performance  in  low  ones.  The 
tracking  community  is  making  progress  toward  seamless  detection  and  tracking 
and  recovering  some  of  this  lost  performance. 

In  this  talk,  we  present  a  method  called  likelihood  ratio  detection  and  tracking 
(LRDT)  that  is  a  step  toward  providing  an  integrated  approach  to  detection  and 
tracking.  We  provide  examples  of  the  application  of  LRDT  to  sonar  and  radar 
detection  and  tracking.  LRDT  is  a  recursive  Bayesian  version  of  track-before- 
detect.  In  LRDT  one  specifies  a  surveillance  region  that  has  a  prior  probability 
less  than  one  of  containing  a  target.  There  is  a  probabilistic  motion  model  that 
specifies  target  motion  within  the  region  as  well  as  the  possibility  of  transiting 
into  and  out  of  the  region.  As  sensor  information  is  received,  it  is  converted  into 
a  measurement  likelihood  ratio  function  and  combined  with  the  prior  likelihood 
ratio  surface  to  produce  a  posterior  surface.  Peaks  in  this  surface  are  used  to 
determine  whether  a  target  is  present  and  to  provide  an  estimate  of  its  state 
(track).  The  process  is  recursive  with  the  posterior  surface  from  one  time  period 
being  updated  for  target  motion  to  become  the  prior  for  the  next  measurement 
likelihood  ratio  function. 
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Outline  of  Talk 


■  Tracking 

□  Bayesian  Tracking 

■  Bayesian  Detection  and  Tracking 

Likelihood  Ratio  Detection  and  Tracking 

■  Example 

□  Passive  Acoustic  Array 

■  Discussion 


Bayesian  Tracking  Framework 


Prior  Distribution 

Start  with  a  probability  distribution  on  the  state  of  the  target 

Account  for  Target  Motion 

□  Specify  distributions  on  motion  and  state  transitions. 

□  Usually  a  Markov  process 

Likelihood  Functions 

Convert  observations  (contacts)  into  likelihood  functions 

Posterior  Distribution 

Combine  likelihood  with  prior  to  compute  the  posterior  distribution 
the  target  state. 

The  posterior  distribution  is  the  output  of  a  Bayesian  tracker 


Likelihood  Functions 


■  The  likelihood  function  converts  an  observation  to  a  function  on 
the  target  state  space: 

L(jl1y)  =  Pr{7  =  y\X(t)  =  s}  for seS 


■  The  observation  (data)  Y=  y  is  known.  The  target  state  5  is  not. 
Thus  L(yU)  is  a  function  on  S,  the  target  state  space 

■  Usually  L(yL)  is  not  a  probability  distribution  on  S 

■  Likelihood  functions  are  also  used  for  association  and 
classification 
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Example 

Bearing  and  Detection  Observations 


Bearing  Likelihood 


Detection  Likelihood 


5 


Bayes  Markov  Single  Target  Recursion 


Recursion 


Initial  Distribution:  p(t0,s0)  for  s0  e  S  at  time  t0 
For  k  >  1  and  sk  G  S 


Motion  Update: 
Calc  Likelihood: 


p  (4  >  st )  =  Js  ft  (ft  I  ft-i )  p(4-i  .  ft-i )  dsh~i 

L{yk\sk)  =  Vv{Yk=yk\X[tk)  =  sk} 


Information  Update: 


P(tk,sh)  =  -^L(Vk\sk)P  (4 > ft ) 


■  Kalman  Filtering  is  a  special  case  of  this  recursion  that  holds  under 
linear  Gaussian  assumptions 


6 


Likelihood  Ratio  Detection  and  Tracking  (LRT) 
vs  Tracking 


■  Tracking 

□  Assumes  target  is  present 

Uses  only  sensor  responses  that  are  above  threshold  (contacts) 

□  Uses  these  responses  to  estimate  state  of  target 

■  Likelihood  Ratio  Detection  and  Tracking  (LRT) 

□  Does  not  assume  target  present 

□  Uses  below  threshold  sensor  responses 

□  Determines 

■  Whether  target  present 

■  Target  state  if  present 

□  Is  a  Bayesian  form  of  Track-Before-Detect  (TBD)  processing 
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Mathematical  Formalism  for  LRT 


■  Same  as  Bayesian  tracking  except 

We  extend  the  state  space  S  by  adding  the  null  state  $  to  represent  the 
possibility  that  no  target  is  present  in  the  area  of  interest. 

■  We  let  S+=Su<j)  be  this  extended  state  space. 

■  We  assume  there  is  at  most  one  target  in  the  region  so  that 

p(  0,  </>)  +  /  p(°>  s)ds  =  1 

s 

■  We  define  the  cumulative  likelihood  ratio  as 


A(t,s) 


P(t,  s ) 

p(t,  <(>) 


Pr{X(i)  =  s  |  Observations  to  time  t} 
Pr{X(J)  =  cj>  |  Observations  to  time  t} 
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Mathematical  Formalism  Continued 


■  Measurement  likelihood  ratio  for  the  observation  Yk  =  y 

r  U,  I  ^  LMA  Pr{yt=y|X(Q  =  4 
^[yl>  Lk(y  10)  Pr  \Yt  =  y  \  X(t)  =  </>} 

■  This  is  the  ratio  of  the  likelihood  of  obtaining  the  observation  Yk  =  y 
given  target  present  at  s  to  the  likelihood  of  obtaining  the  observation 
given  no  target  present. 
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Simplified  Likelihood  Ratio  Recursion: 

Probability  mass  moving  into  region  equals  amount  moving  out 

(Initial  likelihood  ratio)  A(0,s)=  for  sgS 


For  k  >  1  and  s  eS, 

(Motion  Update) 

A  (tt ,  s)  =  qt  (s  1  </>)  +  J  q(s  lit_,  )A (/*_, ,  )  dsk_, 

s 

(Information  Update) 

C(y  is  )=L^y^ 

A t(tt  \sk)  =  C(yk  \st)A-(tt,s) 

(Logarithm  Form) 

In  A  k(tk  lst)  =  ln£(3't  Is*)  +  ln  A~(tk,s) 
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LRT  Implementation  Schematic 


Target_ 
T  rajectory 


Prior  Likelihood  Ratio 


Motion  Update 


Motion  Updated 
Likelihood  Ratio 


T ki  ir\  a  ^ 


Measurement  Likelihood  Ratio  Posterior  Likelihood  Ratio 

Information  Update 
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Velocity  Sheet  Example 


Kinematic  State  Space:  V/elocity  Sheets 
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Passive  Acoustic  Towed  Array  Example 

Simulated  Data- Idealized  situation 

■  Twin  line  array  with  50  phones  in  each  line 

□  75  Hz  design  frequency 

■  Targets  radiate  at  5  frequencies:  50,  55,  60,  65,  and  70  Hz 

■  Background  noise  is  complex  Gaussian  at  phones 

□  SNR  =  -  9  dB  at  each  frequency  at  the  phones 

■  Range  focused  beam  forming 

720  bearings  and  7  ranges:  1000,  1750,  2500,  3250,  4000,  4750,  and  5500m 

■  Measurements  are  mod  square  of  output  from  each  beam  and 
frequency  every  4  seconds 
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Measurement  Log  Likelihood  Ratio  Function 


H0:  Pr{|s(^,w)|2  =  y}  = 


1 


<r„(Ku) 


H{.  yr[\B(k,co)\2  =  j}  = 


exp 


1 


■y 


Kcrn(k,bj) 


exp 


■y 


<j2(k,u)  +  a2(k,uj)  a2  (k,  w)  +  a2n  ( k,u ) 

Noise  a2n  estimated  for  each  beam,  frequency,  and  time 
Signal  a3.  =  .  125cr^  (Alternative:  specify  prob  dist  on  signal  level) 


\njC(\B(k,u)f  I  k 


|  B(k,uj) 

2  |  B(k,u) 

2 

1  n 

cr2(k,uj)  +  cr2n(k,uj) 

a2(k,uf)  a2s(k,u>)  +  al 

,  111 

:(*,«) 

,  °;,(k,uj) 
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Time 


Beam  Output:  Summed  over  Frequency  -  Range  5500m 

Beam  Power  (dB) 
range  5500m,  frequency  50-70Hz 
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Time 


Measurement  Log  Likelihood  Ratio  -  Range  5500m 
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position  north  (km) 


Log  Likelihood  Ratio  Surfaces  in  Geographic  Space 


Motion  Updated  LLR  Surface 
00:00:20 


Measurement  LLR  Surface 
00:00:20 
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Measurement  LLR  at  00:00:20 
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Cumulative  LLR  Marginal  at  00:00:20 


position  north  (km) 


Cumulative  Log  Likelihood  Ratio  Surfaces 


Before  Carve  Out  at  00:00:20 


After  Carve  Out  at  00:00:20 


Cumulative  LLR  Marginal  at  00:00:20 
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Creation  of  Nodestar  Tracks 
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Nodestar  Tracks 


Nodestar  Marginals  at  00:00:00 


Nodestar  Bearing  Tracks 

NS  Tracks  on  Measurement  LLR 
range  5500m 
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Discussion 


■  LRT  example  illustrates  process  of  going  from  unthresholded  sensor 
output  to  track  initiation  without  seams 

■  Process  of  “carving  out”  contacts  and  sending  to  a  tracker  is  a  seam 
which  needs  to  be  addressed  to  complete  seamless  process 

■  The  basic  methodology  shown  in  the  example  has  been  applied  to 
data  collected  at  sea. 

□  System  (LRT  plus  Nodestar)  was  able  to  track  a  large  number  (40+)  of 
targets  simultaneously  and  automatically. 

Had  some  mis-association  problems  when  tracks  crossed  in  bearing 
■  Adding  spectrum  to  target  state  may  improve  this 

□  Has  been  applied  to  matched  field  data  taken  at  sea. 

Good  array  element  localization  and  estimates  of  background  noise  crucial 

■  Further  discussion  of  Likelihood  Ratio  Detection  and  Tracking  may  be 
found  in  Bayesian  Multiple  Target  Tracking  by  Stone,  Barlow,  and 
Corwin,  Artech  1999 
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